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ABSTRACT 

We determine the mass distribution of a spiral galaxy, 2237+0305 using both 
gravitational lensing and dynamical constraints. We find that lensing can break the 
disc-halo degeneracy. 2237-1-0305 has a sub-maximal disc, contributing 57±3 per cent 
of the rotational support at the disc maximum. The disc mass-to-light ratio is 1.1±0.2 
in the I-band and the bulge, 2.9±0.5. The dark matter halo, modelled as a softened 
isothermal sphere, has a large core radius (13.4±0.4 kpc = lArd) to high accuracy 
for the best-fit solution. The image positions are reasonably well fitted, but require 
further rotation information to obtain a unique solution. 

Key words: dark matter - gravitational lensing - galactic dynamics. 



1 INTRODUCTION 

There have been countless papers in recent years discussing 
the matter distribution in galaxies and, in particular, the 
contributions of both the stellar disc and the dark matter 
halo. Any attempt to model the structure of a galaxy with 
multiple components is limited by the information provided 
by the light we receive. This introduces degeneracies in the 
mass distribution when we compare the few constrained pa- 
rameters with the large number of unknowns. 

One technique that has been used to side-step this prob- 
lem is to assume a maximal disc (quantified in Sackett 1997) , 
whereby the disc contributes the majority of the rotation 
(75-95%) at the radius of its maximum circular speed. This 
definition takes into account the contribution to the inner 
rotation curve of a bulge component. Various studies have 
challenged and s upp orted this view and they will be dis- 



cussed in Section 5.3 



Observations of rotation curves of galaxies gave the first 
hints of the existence of dark matter. In the 1960s and 70s, 
extensive studies of the circular rotation of spiral galaxies 
suggested the need for dark matter to account for the miss- 
ing dynamical mass (e.g. see Rubin et al. 1962; Rubin, Ford 
& Thonnard 1978; Faber & GaUagher 1979). The fiatness 
of rotation curves beyond the optical edge of these spirals 
contradicted the expected Keplerian fall-off. In the 1980s 
and 90s, many groups used large scale structure to deduce 
the expected properties of dark matter particles and applied 
these to N-body simulations (Peebles 1984; Navarro, Frenk 
& White 1996). The subsequent calculated profiles can be 
compared with the results of observations for consistency. 

The firm determination of a general trend in the fea- 
tures of rotation curves, as distinct from a few abnormal 
systems, requires the accumulation of data from many galax- 



ies with varying morphological types. Salucci & Burkert 
(2000) constructed 'universal rotation curves' (URCs) with 
luminosity as the only free parameter, from observations of 
~1100 rotation curves. From these URCs, they were able 
to subtract constant mass-to-light ratio discs from the sur- 
face brightness distributions of observed galaxies, and derive 
the expected contribution from dark matter. Salucci (2001) 
found the dark matter haloes were required to have large 
core regions in order to fit the universal curves. Core radii 
of 3-4 disc scale lengths were found to be consistent with 
observations for large spiral galaxies. 

The use of gravitational lensing to disassemble a galaxy 
has been used recently by Mailer et al. (2000). Their anal- 
ysis of the doubly imaged system B 1600+434 did not have 
a unique solution due to the necessity of invoking the TuUy 
Fisher relation for rotational information, the large errors 
on the position angle of the major components and the per- 
turbing effect of a nearby galaxy. In addition, the constraint 
of only two images limits the number of known parameters. 
Their work, however, illustrated the utilisation of lensing to 
break the disc/halo degeneracy. They found a high proba- 
bility for a sub-maximal disc, and the need for some degree 
of a constant density core in the centre of the dark matter 
halo. These results pave the way for our work and demon- 
strate the ability to solve for the mass distribution given a 
reasonable number of known parameters. 

Galaxy 2237+0305 is improbably close with the source 
almost exactly aligned with the macrolensing galaxy. Its 
four, almost symmetric images of the background quasar 
nicely straddle the central nucleus, placing constraints on 
the structure of the inner-most regions of the galaxy. Tradi- 
tionally, it has been difficult to obtain information in these 
regions due to the degeneracy introduced by multiple mass 
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components and smearing of spectroscopic data from non- 
circular stellar and gas motions. Lying at low redshift, it 
displays a visibly extended spiral disc, allowing a detailed 
investigation of its mass distribution beyond the lensing re- 
gions. In addition, two points on the rotation curve have 
been obtained from neutral hydrogen observations (Barnes 
et al. 1999). The combination of its proximity and its im- 
probable alignment with the background quasar makes this 
system the ideal laboratory for studying galactic structure. 

Section 2 of this paper discusses the basics of lensing. 
Section 3 introduces the galaxy 2237-1-0305 and Section 4 the 
mass models, our method and the constraints used. Section 
5 presents our results and their implications for the structure 
of the galaxy. 



2 LENSING 

Gravitational lensing can be treated as a purely geometric 
consequence of General Relativity. If a background source is 
sufficiently close to the optic axis linking the lens and the ob- 
server (approximately within one Einstein radius) , the paths 
of the light rays received from it will have been altered by 
the mass of the lensing object. If the simplest case of point 
masses for both lens and source is assumed, the positions of 
the images observed in the lens plane can be obtained an- 
alytically using simple geometry. This geometry is directly 
generalisable for more complex mass distributions, where ex- 
tended lenses comprise multiple components. The lens equa- 
tion linking these can be written easily. 



where 



(1) 



(2) 



is the reduced bending angle, 6 is the image position, (3 is 
the source position, and all distances are angular diameter 
distances (Narayan & Bartelmann 1999). Throughout, we 
have assumed an Einstein-de Sitter universe with a Hub- 
ble constant of Ho — 70 kms~^Mpc~^, although cosmology 
makes little difference to physical scales of the galaxy due 
to its proximity to the observer. 

The angle of deflection of a light ray travelling from 
source to observer is dependent upon the projected surface 
mass distribution in the lens plane of the lensing object. 
This incorporates both the convergence (the greater the sur- 
face mass density, the greater the deflection) and shear (the 
greater the anisotropy of the lens, the greater the deviation 
of the images from a uniform, circular conflguration) of the 
lens. The bending angle, a is given by. 
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where E is the surface mass density and the integral extends 
over the surface of the distribution (primed co-ordinates). 
The addition of multiple mass components is trivial with this 
formalism and allows a complicated system to be applied to 
lensing. 



Bar Axis (39 deg) 




Image Axis (67 deg) 




^ O 


B 

o 


Galactic Axis (77 deg) Galactic Centre 




• 


dO ^^^^ 


o 

Image A 





Figure 1. Locations of the four lensed images as measured by 
Crane et al. (1991) and the major axes of the principal galactic 
components as measured by Yee (1988). 



3 2237+0305 

2237+0305, a barred spiral Sab type galaxy at a redshift of 
2=0.0394 was first discovered by Huchra et al. (1985), but 
the images were first seen by Yee (1988). Almost collinear 
with the centre of the galaxy (the optic axis) is a background 
quasar at z=1.695. The geometric conditions of this system 
are favourable to macrolensing of the quasar through the 
bulge of the galaxy. 

Huchra's lens, as it is often termed, has been exten- 
sively studied over the past fifteen years, mainly due to its 
proximity to us, and the ability to image the galaxy with 
ground-based telescopes. For this reason, it is the ideal lab- 
oratory in which to study lensing, and more importantly for 
this work, the structure of the lensing galaxy. The lens has a 
visible central bulge, a stellar disc and a perturbing bar. The 
four images straddle the nucleus in an Einstein Cross config- 
uration and their positions have been measured accurately 
(Crane et al. 1991). The major axis of the disc and bulge 
have been measured by Yee at a position angle of 77° , and 
the images at an angle of 67°. The rotation of the images 
away from the major axis of the galaxy is due to the torque- 
ing effect of the bar at a position angle of 39°. Schmidt et 
al. (1998) extensively modelled the mass distribution in the 
bar and its effect on the positions of the images. The an- 
gular positions of the major components of the galaxy are 
displayed in Figure |l| 

The scale lengths of the luminous parts of the galaxy 
have been measured from HST light profiles by Schmidt 
(1996), and their ellipticities by Racine (1991) and Irwin 
et al. (1989), and thus many of the galactic parameters have 
been previously measured. 

Qualitatively, the majority of the convergence will be 
provided by the central bulge (in the absence of a cuspy 
dark matter halo) as it dominates in the region in which 
the lensing is most affected. The disc will contribute both 



© 0000 RAS, MNRAS 000, 



Dissecting a galaxy: mass distribution of 2237+0305 3 



globally and add significant convergence to the inner regions. 
The bar will add little convergence for the lensing, but will 
provide shear through its rotated position angle. In addition, 
the inclination of the disc and flattening of the bulge will 
also contribute to the twisting of the image positions. The 
presence of a dark matter component to the galaxy will add 
convergence and, if non-spherical, shear. Given the highly 
elliptical nature of the visible components already discussed, 
we will assume a spherical halo. 



4 MODELS 

4.1 Mass Models 

We will construct models for the mass distributions of the 
bulge, disc, bar and dark matter halo using parameters from 
previous studies, and map light rays from the measured 
image positions to the source plane. Varying the contribu- 
tion of each component will vary the convergence and shear 
within the images and shift their back-mapped positions in 
the source plane. A potential solution is obtained when a 
particular addition of the four components produces a com- 
mon source position. Clearly, the four images originate from 
the same point in the source plane. The actual position of 
the source quasar cannot be constrained. The centre of the 
galaxy will be considered fixed given the relatively few de- 
grees of freedom. 

Rotation curves for potential solutions will be produced 
and compared with neutral hydrogen rotation measurements 
and the measured mass lying within the images. The combi- 
nation of both lensing and dynamical constraints increases 
the number of constrained parameters and consequently re- 
duces the number of unknowns. 

The models used for the mass distributions of the four 
principal galactic components are standard profiles from the 
literature, tailored to suit this galaxy. 

The bulge is modelled as both a modified de Vau- 
couleurs surface mass distribution (de Vaucouleurs 1948, 
1959), where it is assumed the mass follows the light (con- 
stant mass-to-light ratio) and an exponential surface mass 
profile, as in the models of Schmidt (1996). The modifica- 
tion to both allows the introduction of an ellipticity, e, such 
that for the de Vaucouleurs profile. 



log E(a::, y) = log Eo - 3.33 



Tb 



(4) 



where E is the value of the surface mass density at that 
position, rt is the characteristic scale length of the bulge, 
and e is defined by, 



e= 1- -, 
a 



(5) 



where a and b are the semi-major and minor axes respec- 
tively. The central surface density (EqxIO'^'^^) is denoted 
'bg'. The exponential profile is modelled simply by intro- 
ducing the ellipticity (assumed to be a projection effect). 



T.{x,y) = bge 



(l-e)2 



(6) 
(7) 



and again 'bg' is the central surface mass density. 

The disc is modelled as an exponential surface den- 
sity. Unlike the bulge which is treated with the ellipticity as 
measured, the disc is rotated to its measured inclination of 
1=60° . This involves projecting the volume to a surface mass 
density by rotating the a-axis and redefining co-ordinates. 
If we assume the disc is uniformly distributed in the z di- 
rection, then we can simply write. 



(8) 



where the proportionality includes a factor reflecting the 
thickness of the disc, assumed to be constant, and ra is the 
characteristic disc scale length. Upon rotation about the x- 
axis (such that it becomes the major axis of the ellipse) by 
the inclination angle, i=60°, the surface mass density is the 
integral through the rotated axis, z' , 



S(a;',y') = dc / e 



where the limits of integration bound the original constant 
disc thickness at the inclination angle (taken as A2;=500pc), 
the primed co-ordinates represent the new, observed Carte- 
sian system and 'dc' denotes the central surface mass den- 
sity. 

The bar has been extensively modelled by Schmidt 
(1996) and we will use his surface mass distribution and 
position angle. Schmidt uses a Ferrers model with an ellip- 
ticity, e, 

A 

(10) 



E(.,y) = br(l-^!-|! 



where 'br' is the central surface density, A is the Ferrers 
exponent, and the {x, y) co-ordinates lie in the rotated frame 
of the bar. Schmidt finds different exponents, ellipticities 
and scale lengths depending on the profiles used to fit to 
the light distribution. For an exponential bulge and disc, he 
finds A=2, e=0.64 and b=1.0±0.3 arcsec fit the observations 
best. For a de Vaucouleurs bulge and exponential disc, he 
finds A=0.5, e=0.89 and b=3.1±0.9 arcsec. In our analysis, 
the central surface density (essentially the M/L) will remain 
a free parameter. Here Ecr is the critical surface density 
defined by the geometry, 

_ c' Ds 

The profile of the dark matter halo is controversial. Con- 
ventional theories where ACDM is the preferred cosmology 
have been very successful in explaining the observed large 
scale structure of the universe (eg Peebles 1984). Navarro, 
Frenk and White (1996; hereafter NFW) used N-body simu- 
lations to derive a density profile for such a cosmology. One 
feature of the NFW profile is a cuspy central region with a ~ 
-1 where p ~ r*^. This steepens to a ~ -3 for r^rh, where 
is the characteristic scale length of the halo. Recently, more 
accurate simulations have pushed the cuspiness at r ~ to 
a ~ -1.5 (eg Moore et al. 1999a). Recent observational work 
by de Blok et al. (2001) on low surface brightness galaxies 
has shown the need for a core in the dark matter haloes, us- 
ing optical rotation curves and fitting minimal discs. CDM 
has been challenged further by the missing satellite (Moore 
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et al. 1999b) and angular momentum problems (see recent 
discussion by Sommer-Larsen & Dolgov 2001). The former 
considers the few satellite galaxies observed in orbit around 
our Galaxy compared with that predicted by the theory, 
and the latter refers to the predictions of CDM of too much 
angular momentum loss to support observed disc galaxies. 

Since there is no general agreement we have chosen to 
model the dark matter halo as a softened isothermal sphere 
(Kormann, Schneider & Bartelmann, 1994), 



where av is the velocity dispersion and rc is the core or break 
radius. There are several reasons for this choice of profile. 
Firstly, it has one less parameter than an NFW+core pro- 
file (no scale length as distinct from the core length). Given 
the small number of degrees of freedom we have, an extra 
one is useful. Secondly, the profile naturally asymptotes to 
a flat rotation curve for r rc, and finally, recent results 
find that the two models are indistinguishable when the re- 
maining mass components are taken into account (Weiner, 
Sellwood & Wilhams, 2001, modelled NGC 4123 with both 
profiles and found their shape had a minor effect on the 
results). In addition to this, the choice of a spherical halo 
is a simplification in that it precludes the need to intro- 
duce another parameter. This choice is partially justified by 
the ellipticities in the other components. They are adequate 
to produce the required ellipticity. Furthermore, recent re- 
sults from Ibata et al. (2001) suggest dark matter haloes 
are spherical. They used evidence of the tidal stream from 
the Sagittarius dwarf galaxy to show that the halo potential 
cannot be flatter than q < 0.7 and probably has q > 0.9, 
where q is the axis ratio. 

4.2 Constraints and Previous Work 

High resolution imaging and detailed modelling of Huchra's 
Lens has provided us with many of the model parameters re- 
quired for fitting a mass distribution. Schmidt (1996) used 
I-band HST imaging to measure the disc and bulge scale 
lengths, as well as to determine their ellipticities. These 
results were updated values from previous studies by Yee 
(1988) and Huchra et al. (1985). These data leave the disc 
and bulge mass-to-light ratios as the only fitting parameters 
for these components. The bar has a measured ellipticity, 
position angle and major axis, and within mass-to-light ra- 
tio uncertainties is completely determined by Schmidt. 

Schmidt deconvolved the HST light distribution using 
two different bulge models, an exponential and a de Vau- 
couleurs surface mass profile. Results from the literature 
show that one of these profiles commonly fits a bulge well. 
CaroUo et al. (2001) show from HST observations that the 
bulge mass profile is mildly morphologically dependent. For 
Sa-Sb galaxies, such as 2237-1-0305, both types are observed 
and thus an exploration of both models is prudent. 

One further constraint and one check will be applied to 
the final rotation curve. Neutral hydrogen observations by 
Barnes et al. (1999) at the VLA provide two rotation points 
in the outer regions of the galaxy. In these regions, the visi- 
ble matter has fallen below observational levels and the HI 
is acting as a tracer of the dark matter distribution. Unfor- 
tunately, the data are not of high enough angular resolution 



to probe the rotation in the inner regions of the galaxy. An 
additional piece of information is provided by gravitational 
lensing, where the position of the images in combination 
with the geometry of the source-lens-observer system gives 
the projected mass enclosed within the images (Rix, Schnei- 
der & Bahcall, 1992; Wambsganss & Paczynski, 1994). The 
consistency of this value with the mass distribution found 
will act as a check on the result. We define the image radius 
as the average of the radii of the four images from the cen- 
tre of the galaxy. This corresponds to rim = 0.9 arcsec = 
QlOh^g parsecs. The halo is completely unconstrained ob- 
servationally. We will fit to both the scale length and the 
halo normalisation. 

We therefore have seven parameters for which a fit is 
required - four mass-to-light ratios, one core radius and two 
co-ordinates of the source position. We have as constraints, 
eight co-ordinates of image positions, two rotation points 
and one mass enclosed within the images. These observa- 
tional constraints are given in Table ^. Thus, we have four 
degrees of freedom. 

The galaxy has been previously modelled by many 
groups. Huchra et al. (1985) undertook the initial work on 
the system, measuring ellipticities and scale lengths and pro- 
viding a rudimentary lensing analysis. They assumed a cir- 
cular lens model and showed that the inferred mass-to-light 
ratio is within current values for nearby galaxies. Kent & 
Falco (1988) approximated the galaxy as an oblate spheroid, 
citing the bulge and the bar as the two primary lensing com- 
ponents. Their analysis attempted to fit the observed quasar 
fiuxes to their model and was reasonably successful. Contem- 
poraneously, Schneider et al. (1988) used a single, elliptical 
de Vaucouleurs bulge to model the galaxy mass distribution. 
They fitted two free parameters, the mass-to-light ratio (as- 
sumed to be constant) and the source position, given the im- 
age positions, by minimising the deviation of predicted po- 
sitions to those observed. Unfortunately, they were limited 
by imaging taken with large PSFs and by the simplicity of 
their modelling. Both Kent & Falco, and Schneider's group 
find different positions for the background source, as does 
Schmidt (1996) in his bar-centric analysis. Schmidt com- 
bines the disc and bulge into a single component with the 
ellipticity Racine (1991) measured for the bulge alone. The 
bar was then studied in detail using high resolution HST 
(WF/PC-1 ) I-band imaging and the source position again 
back-mapped from the measured image positions. This is the 
first study to include the dark matter halo in the modelling. 

In order to have a reasonable chance of understanding 
the complex structure of this lens, and hence its lensing char- 
acteristics, one needs to model each component of the lens 
carefully. With the accumulation of data on this galaxy over 
the past fifteen years, this is possible. 



5 RESULTS 

5.1 Source Position and Rotation Curve 

The lensing galaxy was constructed as a two-dimensional 
array of surface mass elements, with an exponential spatial 
scale. This feature allowed the important central regions to 
be highly resolved, while keeping the array dimensions com- 
putationally tractable. The galaxy was rotated from its po- 
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Quantity 



Value 



Source 



Halo 



No constraints 



Bulge 



Position angle 


77° 


Yee (1988) 


EUipticity 


0.31 


Racine (1991) 


dV 






Scale length 


4.1±0.4" (S.lkpc) 


Schmidt (1996) 


Exp 






Scale length 


0.59±0.03" (0.45kpc) 


Schmidt (1996) 


Disc 


Position angle 


77° 


Yee (1988) 


Inclin. angle 


60° 


Irwin et al. (1989) 


dV 






Scale length 


11.3±1.2" (8.6kpc) 


Schmidt (1996) 


Exp 






Scale length 


5.6±0.4" (4.27kpc) 


Schmidt (1996) 


Bar 


Position angle 


39° 


Yee (1988) 


dV 






Sfc. Brightness 


/o=20.4±0.2 


Schmidt (1996) 


EUipticity 


0.64 




Exp 






Sfc. Brightness 


/o=19.9±0.2 


Schmidt (1996) 


EUipticity 


0.89 





Images 



Mass Enclosed 
(xIQIOMq) 


(1.48±0.01) h-^ 
= (1.59±0.01) 


WP (1994) 


Positions (") 


ARA ADec. 


Crane et al. (1991) 


Image A 


0.093 -0.936 




Image B 


0.579 0.737 




Image C 


-0.719 0.266 




Image D 


0.761 -0.419 





Rotation:Vc 



22 ± 1 kpc 
29 ± 1 kpc 



310±15 kms" 
295±15 kms" 



Barnes et al. (1999) 



Table 1. Assumed values for the mass distributions from pre- 
vious work. Image positions are relative to the galactic centre, 
and have positional uncertainties of ±5mas. WP (1994) denotes 
Wambsganss & Paczynski (1994), and the central I-band magni- 
tudes (/q) ''■re measured in mag./arcsec^. dV and Exp denote 
models with a de Vaucouleurs and an exponential bulge, respec- 
tively. 



sition angle to make the x-axis coincident with the galac- 
tic major axis, and the image positions rotated to a right- 
handed co-ordinate system. The four mass components were 
constructed and overlayed on the mass array, forming a sur- 
face mass distribution for the galaxy. Alterations were made 
to this array when varying the halo core radius, Vc or the 
scale lengths of the bulge and disc. The bending angles for 
the particular distribution were then calculated. 

We have calculated bending angles for particular values 
of bulge and disc scale length, for a given bulge profile. These 
values for the four components then were added in linear 
combination varying the mass-to-light ratios of each compo- 



nent. Since the bending angles scale linearly with mass-to- 
light ratio, variations in the overall scaling parameters can 
be considered after the angles have been calculated. 

Rotation curves were carefully constructed to incorpo- 
rate the lack of spherical symmetry in the disc and bulge. 
The functional form for the disc rotation was obtained for 
an exponential disc from Binney & Tremaine (1987). The 
de Vaucouleurs bulge was de-projected into an elliptical vol- 
ume mass density using the profile of Fillmore (1986) and 
the rotation calculated using mass enclosed within a radius 
and equating the centripetal and gravitational accelerations. 

The mass measured by Wambsganss & Paczynski 
(1994) within the images uses the relations of lensing and 
determines the projected mass in a cylinder bounded by the 
images and integrated along the line-of-sight. A suitable dis- 
tribution of mass profiles must satisfy this constraint. The 
consistency between this mass and a solution that repro- 
duces the correct image positions does not provide a com- 
pletely independent constraint as it has been derived from 
lensing. It only demonstrates the true convergence and shear 
that has already been accounted for in the source finding 
program. As such, the consistency of these points provides a 
check that the source finding program is reconstructing the 
correct shear and convergence. 

The process of deciding which distributions adequately 
fit all of the constraints reduces to the minimisation of the 
statistic. For a particular set of profiles for the luminous 
mass components (as given in Table a search through five 
dimensional parameter space (four mass-to-light ratios and 
one core radius) was performed. For each model providing 
a suitable mass enclosed within the images, the eight image 
positions and the two points at the pertinent positions on the 
rotation curve, were calculated. These were used to calculate 
the statistic and the minimum found. In addition to the 
variation of the mass profiles (de Vaucouleurs (dV) versus 
exponential (Exp) bulge, and subsequent changes in disc and 
bar profiles), the scale lengths of the disc and bulge were 
varied to their la observational limits in order to obtain a 
more complete representation of the parameter space. 

The incomplete Gamma Function provides the goodness 
of fit for a system with n degrees of freedom (d.o.f.) and a 
statistic. For a system with four d.o.f., a solution with 
the probability to la (68%) of being real (that is, not by 
chance) is given by solving. 



F(n/2;x72) =0.1 



(13) 



where F{a; (i) is the incomplete gamma function (Numer- 
ical Recipes). The solution for n=4 is Xmax=2.3. This is 
a strict value from the formal theory and requires a high 
degree of accuracy to be attained. Considering the idealistic 
mass models we are using and the small number of d.o.f., we 
would not expect to reach this low value for the x^ . Instead, 
the minimum value found will be taken as an adequate fit, 
provided it does not exceed this value excessively. The error 
analysis for the best fit result will provide information on 
the sensitivity of the solution on the different parameters. 
We will take the errors on a parameter to be where the x^ 
increases by one while holding all other parameters constant. 

The results for the model with a de Vaucouleurs bulge 
is shown in Table |^, and an exponential bulge in Table |^. 

These tables clearly demonstrate the superiority of the 
former over the latter. The x^ parameter for the exponen- 
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n 


3.1 


3.1 


3.1 


2.8 


3.4 


rd 


8.6 


9.5 


7.7 


8.6 


8.6 


Tc 


15.3±0.4 


13.4±0.4 


14.9±0.5 




14.5±0.4 




266±4 


233±4 


246±4 




257±4 


bg 


9.66±0.04 


9.53±0.04 


9.47±0.04 




8.66±0.04 


dc 


242±24 


506±30 


571±38 




447±35 


br 


824±23 


821±20 


818±20 




783±22 




10.95 


5.26 


6.99 


>500 


9.81 



Table 2. Best-fit parameters for a given model, exponential disc 
and de Vaucouleurs bulge. The scale lengths (r;,, r^) are given in 
kpc; is the velocity dispersion (km/s) of the halo profile; bg, 
dc and br are the central densities (Mqpc"'^) for the bulge, disc 
and bar respectively. The bulge values shown are divided by 10^ 
of their actual value. 





0.45 


0.45 


0.45 


0.43 


0.47 


rd 


4.27 


4.58 


3.97 


4.27 


4.27 


Tc 


15.2 


15.2 


15.2 


15.9 


15.7 




243 


243 


243 


262 


239 


hg 


3.0x10'' 


3.0 xlO" 


3.0x10"' 


3.2x10" 


2.9x10"' 


dc 


0.0 


0.0 


0.0 


3.6 


0.0 


br 


924 


924 


924 


934 


894 




140 


140 


140 


106 


210 



Table 3. Same as for Table ^ except for an exponential disc and 
bulge. Symbols have the same meaning, but the bulge scaling is 
correct. Uncertainties were not calculated due to the large 
values. 



tial bulge results is evenly distributed between error in the 
rotation curve points and the image positions. It seems the 
rotation curve cannot be sufficiently produced to satisfy the 
mass enclosed within the images constraint as well as the two 
points measured in the curve. Similarly, the combination of 
bulge, disc and bar cannot produce the required ellipticity 
as well as an acceptable rotation curve. 

Conversely, the results for the de Vaucouleurs bulge and 
exponential disc are quite promising. None of the combi- 
nation of scale lengths produced an acceptable (<2.3), 
however, the values are of the right order of magnitude. The 
best result is produced at the la upper limit of the disc scale 
length as measured by Schmidt. Without a continuity of so- 
lutions between these scale lengths, it is difficult to predict 
if an acceptable solution occurs within the measured limits. 
Without the added accuracy of further rotation points, we 
will take this solution to be acceptable and base our analysis 
upon it. The for this result is distributed evenly between 
error in the rotation points and the image positions. 

Figure ^ shows the rotation curve produced from this 
solution and the components used to construct it. The total 
rotation curve is still rising at the edge of the plot due to 
the influence of the halo. Here, the three luminous compo- 
nents are falling away, particularly the bulge and bar which 
are essentially contributing no mass. This rise could pose 
a problem in that rotation curves are observed to be rea- 
sonably flat or slowly falling. Simply scaling down the halo 
core size would alleviate this problem, however it would also 
affect the convergence at the image positions. Although the 
halo does not contribute significant mass within the images, 
losing this would require a scaling up of another component. 



CoryiponGnt 


Mass Enclosed (x 1(P Mq) 


Halo 


6.8±0.2 


Bulge 


130.3±0.5 


Disc 


7.1±0.4 


Bar 


11.8±0.3 


Total 


156.0±1.4 



Table 4. Projected mass within the image radius divided into 
the four mass components for the best fit solution of Table 



Group 



AN (arcsec) AE (arcsec) 



1. Kent & Falco (1988) -0.02±0.01 

2. Schneider et al. (1988) 0.015±0.005 

3. Schmidt et al. (1998) -0.014 

4. This work 



0.003 
+0.001 
0.014±0.001 



-0.08±0.02 
-0.004±0.005 
-0.063: 



-0.010 
-"4-0.009 

-0.072±0.001 



Table 5. Mean and standard deviation of the source position for 
different groups, as measured from the centre of the galaxy. AN 
and AE refer to the offset from the centre of the galaxy in the 
north and east directions. 



As all of the other modelled components contribute shear, 
this would effect the overall image positions and increase the 
X^- Another halo profile may provide a better fit. 

As mentioned above it is pertinent to calculate the pro- 
jected mass of each of the components within the images. 
This involves the integration of mass within a cylinder out 
to the image radius and through the line-of-sight. Table ^ 
displays the results. The bulge clearly dominates the con- 
vergence within the image region. A percentage variation in 
the overall scalings of the halo, disc and bar components 
will make less difference to the lensing than a change in 
the bulge value. These results demonstrate the nature of 
stacking these mass components. An increase in the halo 
mass (for example, due to a small core region) would neces- 
sitate a decrease in the bulge. For this to be an adequate 
solution, however, the halo cannot be too large, otherwise 
it contributes too much rotation at the HI points, and the 
loss in shear from the bulge must be compensated for in 
another elliptical component. Thus, the combination of re- 
quiring correct convergence and shear is able to break the 
disc/halo degeneracy, removing the need to consider both 
maximal and minimal discs (Mailer et al. 2000). 

In addition to finding the best-fitting bulge and disc 
contributions to the rotation curve, one can perform a sim- 
ilar analysis on the source position for each potential so- 
lution. The results of this analysis are shown in Table ^ 
compared with those from previous groups. The source po- 
sition for a given configuration was found by averaging the 
backmapped positions of the four images. Our results fit 
well with those of Schmidt et al. (1998) and Kent & Falco 
(1988), but not with Schneider et al. (1988). The consistency 
of three of the results is encouraging. 

The image positions reproduced with this model are not 
adequate to definitively conclude this as a solution, however 
they are still extremely close to the observed positions. Fig- 
ure ^ displays the model and real image positions, as well as 
the location of the source. 

The final piece of relevant information that can be de- 
rived from these results is the total mass of each of the galac- 



© 0000 RAS, MNRAS 000, 



Dissecting a galaxy: mass distribution of 2237+0305 7 



400 




Figure 2. Rotation curve for the best solution (ri,=3.1kpc, rrf=9.5kpc, rc=13.4kpc, (T^=233km/s, bg=9.53x 10'Mq/pc^, 
dc=506MQ/pc^, br=821MQ/pc^) as in Table ^. The curves are total (solid), halo (dash-dot), disc (dotted), bulge (short dashed) 
and bar (long dashed). The two vertical lines denote the rotation points from HI measurements and their Ic uncertainties. 



Component 


Mass (<30 kpc) hj^ Mq 


Halo 


(3.7±0.1) X 10" 


Bulge 


(6.7±0.1) X 10"' 


Disc 


(1.7±0.1) X 10" 


Bar 


(2.3±0.1) X 10"' 


Total 


(6.2±0.2) X 10" 



Table 6. Masses of the four galactic components from the poten- 
tial solutions of Table ^. The values given are total mass within a 
radius of 30 kpc in solar masses. The range of halo values demon- 
strate the variation in total halo mass for the potential solutions. 

tic components. The optical disc extends out to r ~ 25-30 
kpc (Yee 1988) and the dark matter halo will presumably 
continue far beyond this radius. As such, the masses calcu- 
lated will be those enclosed within a sphere with radius r 
= 30 kpc, centred on the galactic nucleus. The results are 
displayed in Table ^. These values are consistent in their 
order of magnitude with nominal masses for galactic com- 
ponents. The dark matter halo contributes ~60 per cent of 
the dynamical mass of the galaxy within the optically visible 
region. If the halo was allowed to extend further, as it clearly 
does from its rotation curve, it would contribute much more 
to the overall galactic mass. 



5.2 Mass-to-Light Ratio 

In order to facilitate his analysis of the bar in 22374-0305, 
Robert Schmidt decomposed light profiles by taking cuts 
along the major and minor galactic axes. These profiles can 
be used to investigate the mass-to-light ratio of the bulge 
and disc in the I-band (Schmidt, private communication). A 
deviation from a constant mass-to-light ratio across an indi- 
vidual component could either indicate the non-applicability 
of the mass profiles to this galaxy, or an intrinsic colour vari- 
ation. Such a failure of the hght to trace the mass would have 
implications for analyses of galactic discs where the assump- 
tion of constant mass-to-light ratio is often made. 

The major axis light profile, once calibrated to units of 
solar luminosity, clearly displays the addition of two distinct 
components. The central region is dominated by a steeply 
falling bulge region, which fiattens to a more gradually de- 
creasing disc. The luminosity was calculated assuming a so- 
lar absolute magnitude in the V-band and correcting it to 
the I-band using the Vilnius spectral shifts of Bessell (1990). 

We calculate the mass-to-light ratio for the disc, ini- 
tially, by simply dividing the best fit mass profile by the 
luminosity profile (see Figure ^) . The location of a probable 
spiral arm is visible as a dip at a radius of ~13 kpc. The plot 
is not fiat, as one would expect if the disc hght followed a per- 
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X Offset from Galoctic Centre (orcsec) 



2000 
Radius (pc) 



Figure 3. Comparison of the model and measured image posi- 
tions for the best-fit solution. The crosses denote the model po- 
sitions, the adjacent diamonds the measured, the triangle is the 
assumed galactic centre, and the isolated diamond is the source 
position. The Ict uncertainties on the measured image positions 
correspond approximately to the size of the diamonds. The po- 
sitions are labelled A, C, B, D clockwise from the bottom using 
the Yee (1988) convention. 




.0x10' 1.5x10' 
Radius (pc) 



Figure 4. Disc I-band mass-to-light ratio as found by dividing the 
major-axis luminosity profile into the surface mass. The central 
4 kpc has been removed to eliminate the effect of the bulge. 



feet exponential with the scale length used. Instead, there is 
a noticeable gradient to the ratio. This is possibly an indica- 
tion of further structure beyond the simple models used. An 
additional minor difference between our models and those 
of Schmidt is our use of a thick disc. Finally, the presence 
of the bar affects the underlying mass distribution. Interest- 
ingly, if the disc scale length is increased to ~12kpc, the 
gradient virtually disappears. This does, however, neglect 
the bar influence. As a crude estimate of the value of the 
ratio, however, one can average the result from the curve to 
find. 



Figure 5. I-band mass-to-light ratio profile for the bulge compo- 
nent. The central 700 pc have been omitted to reduce the effect 
of contaminating quasar light. The sharp rise at r ~3.3 kpc is the 
region between the disc and bulge where the mass distributions 
fail to fit the light profile adequately. 



where I indicates the wavelength band and d denotes the 
disc. 

This value is consistent with those from the literature. 
Syer, Mao & Mo (1998) used disc stability arguments to 
place an upper bound on the I-band mass-to-light ratio, 
(M/L)/ < 1.9hioo = l.Shyo. An independent study using 
population synthesis models (Boissier & Prantzos, 2000) 
gives (M/L)/ = 1.0-1.3hioo = 0.7-0.9h7o. Finally, Sommer- 
Larsen & Dolgov (2001) use the I-band TuUy Fisher rela- 
tion of Giovanelli et al. (1997) and their warm dark matter 
simulations to limit the mass-to-light ratio in this band to 
(M/L)/ — 0.6-0.7h6o±io- These independent results fit well 
with the value we find. 

By fitting an exponential curve to the Ught profile, out- 
side of the bulge region of influence, the hght contribution 
of the disc can be removed and the bulge mass-to-light ratio 
calculated. This technique assumes that the disc hght fol- 
lows the exponential decay to the centre of the galaxy and 
is not disrupted by the other mass components in that re- 
gion. We ignore the Ught emitted within the image region, 
as this is contaminated by quasar light. The final profile for 
the best-fitting bulge model (bg = 9.53 x 10^ Mopc"^) is 
displayed in Figure |^. This profile is reasonably flat over a 
large range of radii, but increases sharply for r > 3 kpc. This 
is the beginning of the region where both the disc and bulge 
are infiuential. Here, the disc exponential profile is a poor 
fit to the light profile and the disc subtraction removes too 
much light. This indicates the inadequacy of the two profiles 
to accurately account for the mass. It is possible the inclu- 
sion of the bar would add mass to the inner regions of the 
bulge area (its effectiveness is negligible beyond r ~3 kpc) 
and increase the mass-to-light ratio. This could account for 
the slope at r > 3 kpc. Unfortunately, as the bulge and bar 
occupy the same radial regions, they cannot be deconvolved 
as the disc could. Taking the data from radii, r = 0.7-3.3 
kpc, the average mass-to-light ratio is. 



(M/L)/ 



1.1 ± 0.2/17 



(14) 



(M/L)/,6 = 2.9 ±0.5/170, 



(15) 
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where b denotes the bulge component. Fukugita, Hogan & 
Peebles (1998) find a value of (M/L)i3,6 = e.Silj.o corre- 
sponding to (M/L)/,6 = 5.0±]I[j g in the I-band using synthe- 
sis models and data from other groups. This determination 
just fits with our results. It is encouraging to note the higher 
mass-to-light ratio in the older population bulge stars than 
for the younger disc stars. 



5.3 Maximality 

Maximal discs have been used in many studies to investi- 
gate the dynamics and structure of galaxies, as a means of 
deconvolving the disc from the dark matter halo. Physically, 
maximal discs have been found to be inconsistent with ob- 
servations of gas motions in spiral galaxies combined with 
theoretical models (for example, Kranz, Slyz & Rix 2001). 
Observations of spiral galaxies have challenged the sugges- 
tion that discs are maximal owing to the need for a signif- 
icant dark matter halo (e.g. Bottema 1997). Weiner, Sell- 
wood & Williams (2001) used observed velocity dispersions 
and gas dynamical simulations to reproduce the observed 
bar and spiral arms of a large spiral galaxy. They found the 
disc was maximal (80-100%) to high confidence. Conversely, 
Courteau & Rix (1999) find that their sample HSB galaxies 
are sub-maximal to high confidence using residuals from the 
TuUy-Fisher relation and adiabatic infall of luminous mate- 
rial into dark matter potentials. Bottema (1993) used stellar 
velocity dispersions to infer the contribution from stars to 
the rotation curves in twelve discs. He finds Vstara/vtotai ~ 
63±10 per cent, below maximal according to the definition 
of Sackett (1997). 

The small variation in disc contribution to the 
2237+0305 rotation curve, as demonstrated in Figure]^, pro- 
vides a good determination of the degree of maximality in 
this galaxy. The overall rotation curve, although not con- 
strained observationally in the region where this calculation 
is made (the disc maximum at r ~ 2.2rtj), is also reasonably 
tight given the HI constraints and the profiles we have used. 

The contribution of the disc to the rotation has already 
been determined by its mean central surface mass density 
(dc — 506±30 M0pc~^). This corresponds to a maximum 
rotation oi Vdisc{2.2rd) = 163±5 kms~^. The maximum ro- 
tation is calculated to be vtotai{2.2rd) ~ 288±5 kms~^. The 
percentage contribution of the disc to the rotation, the de- 
gree, is therefore, 

Vdisc{2.2rd) 



Image Flux Ratio (Agol et al.) Flux Ratio (this work) 



57 ± 3% 



(16) 



Vtotai{2.2rd) 

This value fits well with that found by Bottema (1993), 
and is well defined for the potential solutions presented in 
this work. The disc is clearly sub-maximal. 

5.4 Flux Ratios 

The fiux observed from each image in a gravitationally 
lensed system, is a direct measure of the magnification in 
that region of the lens plane. A comparison between the 
fiux ratios of the images observed in 22374-0305 and those 
predicted by the solutions can further act as a check on the 
results. 

Observed fiuxes in individual images are a combina- 
tion of magnification due to microlensing, macrolensing and 



A 
B 
C 
D 



0.9±0.1 

1.0 
0.5±0.1 
0.9±0.1 



0.7±0.1 

1.0 
0.6±0.1 
l.OitO.l 



Table 7. Infrared flux ratios (relative to the B image) of the A, 
C and D lensed images, calculated from the fluxes of Agol et al. 
(2000) and our corresponding results. 



intrinsic variability coupled with time delays. Agol, Jones 
& Blaes (2000), however, have measured IR fiuxes (8.9 & 
11.7/im) for the four components and from these the ratio 
of fiuxes can be calculated. In this region of the spectrum, 
microlensing events are not observed and it is therefore pos- 
tulated that these observations sample an extended region 
of the source. In addition, the infrared fiuxes are not sensi- 
tive to the dust reddening effects of optical light travelling 
through the galaxy. Thus the IR fluxes should measure the 
macro-magnification. These fluxes, relative to the B image, 
are displayed in Table ^. 

The fiux ratios were calculated from our results by tak- 
ing the ratio of the magnification for each of the images rela- 
tive to the B image. The magnification is calculated by tak- 
ing the ratio of areas of triangles around the images mapped 
from the image to the source plane. They were calculated 
from the best-fit solution of Table |2| and are displayed in 
Table |^. Our results are consistent with the observations. 

5.5 Implications for the Dark Matter Halo 

The halo model used is a generic profile that is analytically 
simple but has little physical motivation. It does, however 
provide a mass profile with varying slope, and this is a useful 
attribute if one wishes to study the gradient of the mass 
distribution. 

The best-fit solution for the mass distribution of this 
galaxy is not an adequate fit. The rotation curve is rising 
at the outer regions studied instead of falling, and the im- 
ages are not close enough to their measured points. A halo 
with a smaller core would alleviate the former problem but 
require a change of the other mass distributions to address 
the latter. The mass-to-light profile of the disc suggests an 
increase in the disc scale length. Such a move would change 
the shear introduced by the disc. In any case, the core radii 
found as the best solutions in this study have large core 
radii, Vc ~ 13-16kpc. These values are not consistent with 
the cuspy central regions of the CDM profiles. In order to 
reduce the and find an adequate solution, we firstly re- 
quire further rotation data. The more points we have with 
good accuracy, the more likely we are to be able to elimi- 
nate profiles and find the best solution. Until more data is 
obtained, it is not worth trying to twist the halo profile to fit 
the parameters better. If the error bars on the HI rotation 
points were doubled, the fit would become quite consistent 
with the observations. The Barnes et al. (1999) observations 
were undertaken with the VLA C array, a compact config- 
uration with low angular resolution. More informative data 
is attainable with higher resolution observations. 

The surface mass distribution of the four components 
combined can be represented on a log-log plot to study 
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Figure 6. Log-log plot of the surface mass density of the inner 
regions of the best-fit solution (solid line). The short-dashed line 
indicates the slope of the profile in the central lOpc (S oc r""'^). 
The long-dashed displays the transition to a more isothermal pro- 
file (S oc r~^) outside of the core. 



the slope as a function of radius. This information is dis- 
played in Figure ^ The dashed lines are normalised fits to 
E oc r~°'^ (short-dashed) and E oc r^^ (long-dashed). The 
central slope overall is reasonably steep given the large in- 
fluence of the bulge. The dark matter halo has zero slope in 
the inner regions given its large core radius. The transition 
to isothermal is expected given the dominance of the halo 
in the outer regions. 



6 CONCLUSIONS 

We have undertaken a study of the structure of galaxy 
2237-1-0305 using dynamical and gravitational lensing con- 
straints. The combination of these techniques allows the 
problems related to the disc/halo degeneracy and maximal- 
ity of the disc to be overcome. 

We find the galaxy is adequately modelled with four 
mass components, a central bar, stellar disc, bulge and dark 
matter halo. The bulge contributes ~85 per cent of the con- 
vergence of the lensing, within the image radius, and is thus 
the principal component to which the lensing is sensitive. 
The disc, bulge and bar all contribute shear, while the as- 
sumption of sphericity for the dark matter halo seems ad- 
equate for this analysis with the observational information 
available. 

Within the visible radius of the galaxy (r < SOh^Q^ kpc) , 
the dark matter halo contributes ~60 per cent of the total 
mass. The potential existence of a core region in the softened 
isothermal sphere halo is suggested to high confidence by 
the results. The best reconstruction of the source position 
occurs for a halo with a non-cuspy central region (within 
the image radius), and a core radius of 13.4±0.4kpc, which 
is equivalent to 1.4?-^ for this model. 

The disc is found to be sub-maximal to Scr with an av- 
erage contribution of 57±3 per cent to the rotation at the 
maximum of its circular speed. This is a low result, but 
clearly demonstrates the dominance of the dark matter halo 



in the galactic mass. The disc and bulge I-band mass-to- 
light ratios were calculated using the mean and standard 
deviation of the potential solutions. They were found to be 
(M/L)j,d = l.l±0.2h7o and (M/L)/,^ = 2.9±0.5h7o, respec- 
tively. The flux ratios between the quasar images calculated 
from their magnifications was found to be consistent with 
infrared data. 

With the addition of rotation curve information, the 
dark matter profile can potentially be found uniquely, and 
comparisons with the theories of dark matter particle candi- 
dates undertaken. With a full rotation curve, the dark mat- 
ter profile can be found uniquely. The emission lines Ha, 
[OIII]A5007 and [NII]A6583 are all present in the galaxy 
spectrum (unpublished spectrum from ANU 2.3m, A. Osh- 
lack, private communication) and would be suitable for this 
study. 
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